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On  the  Solutions  of  the  Differential  Equation  yvl  =  xy. 

I .  Analysis* 

R.  L.  Duty**  and  W.  H.  Reid 
Brown  University 
Providence,  Rhode  Island 


Abstract 

The  differential  equation  yvi  =  xy  plays  an  important 
role  in  the  asymptotic  treatment  of  the  stability  of  viscous  flow 
between  contra-rotating  cylinders  and,  in  one  limiting  case, 
solutions  of  this  equation  are  required  that  remain  bounded  as 
x  — >  +oo  .  A  set  of  standard  solutions  have  therefore  been  defined 
such  that  three  of  them,  denoted  by  A^x),  are  bounded  as  x  — >+oo, 
while  the  remaining  three  solutions,  denoted  by  Bj^x),  are 
unbounded  as  x  — >  +oo .  The  contour  integral  representations  of 
these  solutions  are  given,  together  with  their  power-series  and 
asymptotic  expansions.  It  is  also  shown  that  a  slightly  modified 
set  of  these  solutions  provide  a  "numerically  satisfactory"  set 
over  the  entire  interval  -co  <  x  <  +oo  . 
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Office  of  Naval  Research  under  Contract  Nonr  562(07)  with 
Brown  University. 
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1 .  Introduction. 

In  a  recent  study  of  the  stability  of  viscous  flow 
between  rotating  cylinders  (Duty  &  Reid  1961)  it  has  been  shown 
that,  when  the  cylinders  rotate  in  opposite  directions,  the  solu¬ 
tion  can  be  expressed  asymptotically  in  terms  of  the  solutions  of 
the  comparison  equation 

yvl  =  xy.  (1) 

This  equation  clearly  plays  much  the  same  role  in  the  present 
theory  as  Airy's  equation  (y'^xy)  does  in  the  asymptotic  theory 
of  second-order  differential  equations  with  a  simple  turning  point 
and,  once  a  suitable  set  of  solutions  have  been  defined  and 
tabulated,  the  solution  of  the  stability  problem  can  be  made 
explicit.  For  the  present  purposes  it  is  sufficient  to  consider 
only  real  values  of  x  though,  more  generally,  if  allowance  were 
to  be  made  for  the  possibility  of  overstable  modes  of  instability, 
it  would  then  be  necessary  to  extend  the  analysis  to  complex 
values  of  x. 

There  are,  of  course,  many  different  ways  of  defining 
a  standard  set  of  solutions  of  equation  (l) .  Heading  (1957)>  for 
example,  has  recently  studied  a  class  of  differential  equations 
which  includes  equation  (l)  as  a  particular  case.  The  solutions 
defined  by  Heading,  however,  are  not  the  most  convenient  ones  for 
the  present  purposes  since  they  are  not  all  real  for  real  values 
of  x.  We  shall  instead  find  it  convenient  to  define  a  standard 
set  of  solutions  in  such  a  way  that  three  of  them  remain  bounded 
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as  x  — >+aj  (the  remaining  three  being  then  necessarily  unbounded 
as  x  — »  +co ) .  The  solutions  defined  in  this  way  are  thus  analogous 
to  the  solutions  of  Airy's  equations  usually  denoted  by  Ai(x)  and 
Bi(x) . 

The  relevance  of  equation  (l)  to  the  problem  of  the 
stability  of  Couette  flow  was  first  recognized  by  Meksyn  (19^6, 
1961)  who  attempted  am  asymptotic  treatment  of  the  problem  on 
the  hypothesis  that  the  wave-number  of  the  disturbance  (rather 
than  the  Taylor  number)  is  large.  Although  there  are  many  super¬ 
ficial  similarities  between  these  two  approaches,  there  still 
remain  some  important  differences  that  have  not  been  fully 
resolved.  Furthermore,  Meksyn’s  analysis  of  both  the  stability 
problem  and  the  comparison  equation  was  limited  to  the  leading 
terms  in  the  asymptotic  expansions  of  the  solutions  and  so  gave 
no  information  on  their  behavior  near  the  turning  point. 
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2.  Power-series  solutions . 

Six  fundamental  solutions  of  equation  (l)  can  be  found 
In  the  form  of  power-series  by  elementary  methods  (cf .  Heading 
1957 ) .  Ihey  can  conveniently  be  written  In  the  fora 

yn(x)  =  x11"1  2  agx7s  (n=l,2, . . .,6)  (2) 

s=o 

where  a,  =  1 

o 

/  n  \  |  S 

and  as  =  (n+fs-lji  n  (n+7p)  for  s  ^  1. 

P  ® 

We  now  wish  to  determine  those  combinations  of  the  power-series 
solutions  that  remain  bounded  as  x  — >  +00  and  for  this  purpose 
we  turn  to  a  discussion  of  the  solution  of  equation  (l)  by  the 
method  of  contour  Integration. 
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3.  Contour  Integral  solutions. 

The  solution  of  equation  (l)  by  means  of  a  Laplace 
contour  integral  leads  to  the  integral  representation 

y(x)  «  J  exp(xt-t7/7)dt,  (3) 

C 

where  C  is  any  open  contour  that  starts  in  one  and  ends  in  another 
of  the  shaded  sectors  shown  in  figure  1 .  If  we  now  let  Lrs  denote 
a  contour  starting  in  sector  r  and  ending  in  sector  s  then,  since 
the  integrand  of  (3)  has  no  singularity  in  the  finite  t-plane, 
all  contours  for  given  values  of  r  and  s  are  equivalent  by  Cauchy's 
theorem . 

We  Shall  be  particularly  concerned  with  contours  Lrs 
for  which  s=r*-l  and  r=l,2,  ...,7  with  L^g  =  L,^.  For  such  a 
contour  we  will  denote  the  solution  by 

ur(x)  =J  exp(xt-t7/7)dt  (s=rKL).  (4) 

Lrs 

Note  that  in  this  way  we  have  defined  seven  solutions.  But  again, 
by  Cauchy's  theorem,  we  have 

2  u  (x)  *  (£)  exp(xt-t7/7)dt  =  0. 
r=l  r  T 

We  now  wish  to  obtain  the  asymptotic  behavior  of  the 
solutions  (4)  as  x  — >  +oo  for  we  can  then  define  a  set  of  standard 
solutions  that  are  real  when  x  is  real  and  that  have  the  desired 
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4.  The  asymptotic  expansions  of  the  solutions  as  x  — » +oo  . 

In  the  application  of  the  saddle-point  method  (or,  more 
generally,  the  method  of  steepest  descents)  to  the  solutions  (4), 
it  is  convenient  to  let  t  =  so  that 


exp  {x^^w(t)}  dT 


(s=r+l),  (5) 


where 


7 

w(t)  = 


(6) 


The  saddle-points,  which  are  located  at  the  roots  of  w'(t)  =  0, 
therefore  lie  on  the  unit  circle  in  the  T-plane  at  the  positions 
given  by  the  six  roots  of  +1  (see  table  l) .  If  these  saddle- 
points  are  denoted  by  x^  ( 1=0,1, 2, . . .,5)  then  we  may  note  that 

w(Ti)  =  1  x  i  2111(3  w"(t1)  =  -6t*,  (7) 


* 

where 

then 

and 


denotes  the  complex  conjugate  of  If  we  also  let 


w  =  u+iv  and  x  =  5+itj, 

uU,tl)  «  K  -  i  5(£6  -  21? V  +  35*V  -  7ti6) 

v(5,n)  =  ti  +  Y  ^  -  21^2  +  35ti2^4  "  TC6). 


(8) 


(9) 


The  curves  u(5,rj)  =  constant  form  a  family  of  contour  lines  in 
the  x -plane  and  the  curves  v(5,t])  =  constant  are  their  orthogonal 
trajectories  (except  at  the  saddle-points) .  In  particular,  the 
steepest  paths  through  x.  are  given  by 

v(5,n)  ■  v(51,ri1) . 


(10) 
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The  explicit  determination  of  the  contour  lines  and  steepest 
paths  that  pass  through  the  saddle-points  Is  clearly  a  complicated 
matter  and  only  a  qualitative  picture  of  them,  adequate  for  the 
present  purposes,  is  shown  in  figure  2. 

When  a  path  Lrg  (with  s  not  necessarily  equal  to  r+l) 
is  deformed  into  a  path  of  steepest  descent,*  it  may  pass  through 
one  or  more  saddle-points.  The  dominant  contribution  to  the  solu¬ 
tion  associated  with  such  a  path  then  comes  from  the  highest 
saddle-point  on  the  path.  Since  the  elevations  of  the  saddle- 
points  are  determined  by  re{w(xi)},  their  relative  heights  can  be 
ordered  according  to  the  following  scheme: 


t:  +1 
o 

and  +g- 

\  and 


(11) 


Consider  now  the  solution  u,(x)  associated  with  the 
path  L1g .  Biis  path  can  be  directly  deformed  into  a  steepest 
descents  path  passing  through  and,  by  the  saddle-point  method, 
we  then  obtain 


u^x)  -  i(n/3)1//2  x"5/12  exp(-  |  x7/6).  (12) 


It  may  not,  in  fact,  be  so  deformable  and  in  such  cases  an 
equivalent  path  must  be  chosen  that  is  deformable  into  a 
steepest  descents  path. 
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The  descending  series  associated  with  this  solution  could  be 
found,  if  necessary*  by  the  method  of  steepest  descents  but  this 
is  somewhat  difficult  and  not  necessary  for  the  present  purposes. 
Although  we  are  concerned  here  with  real  values  of  x  only,  it 
can  be  shown  (cf .  Heading  1957)  that  the  expansion  (12)  remains 
valid  (in  the  strict  sense)  in  |arg  x|  <  2v/J .  In  the  Poincare 
sense,  however,  it  is  valid  in  the  larger  range  |arg  x|  <  x. 

The  solutions  Ug(x)  and  u^.(x)  associated  with  the  paths 
and  respectively  can  be  treated  in  a  similar  manner  to 
yield 


u2(x)  ~  +  e+lti/3(V3)1/2  x-5/12  exp(+  |  TgX7/6)  (13) 

and  u^(x)  ~  -  e-lll/3( tc/3)1^2  x“^12  exp(+  j  t^x7//^)  .  (l4) 


Since  re^g)  -  re(T^)  =  -  jj,  these  two  solutions  are  dominant 
with  respect  to  u^(x)  but  they  are  subdominant  (though  not 
maximally  so)  with  respect  to  any  solution  whose  path  starts  or 
ends  in  sectors  4,5,  or  6. 

Thus,  from  these  results,  it  is  clear  that  any  set  of 
solutions  that  is  bounded  as  x  — >  +oo  must  be  formed  from  linear 
combinations  of  u^(x),  ug(x),  and  u^(x) .  The  further  requirement 
that  the  solutions  be  real  for  real  values  of  x  then  suggests  the 
choice : 


Vx)  =  2xi  V 

a2(x)  ■  2xi  (u2+u7^ 


> 


(15) 


and 


j 
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The  leading  terms  in  the  asymptotic  expansions  of  these  solu¬ 
tions  are  then  given  by 

Ai(x)  ~  -g-C 3^t ) exp(„  6  x7/6^ 

A2(x)  -  (3 x)-1/2x-5/12  exp(-  |  x7/6)8in(A^2  x7/6  +  |),  ^  (l6) 

and 

A3(x)  -  (aO-VV.5/ 12  exp(-  |  X7/6)  COS  xr/6  + 1) . 


Although  the  choice  (15)  would  appear  to  be  the  most 
natural  one  in  view  of  the  symmetry  that  has  been  achieved, 
these  solutions  do  not,  unfortunately,  form  a  "numerically 
satisfactory"  set  of  solutions  for  large  negative  values  of  x. 
This  matter  will  be  examined  more  fully  in  the  following 
section  where  advantage  will  be  taken  of  the  fact  that  an 
arbitrary  multiple  of  A^(x)  can  be  added  to  Ag(x)  or  A^(x) 
without  changing  the  dominant  series  in  their  asymptotic 
expansions . 

To  complete  the  set  of  standard  solutions  we  must  now 
consider  the  solutions  u^u^u^,  and  Ug.  Consider  first  the 
solution  u^  associated  with  the  path  L^.  This  path  can  be 
deformed  into  a  steepest  descents  path  but  it  is  seen  from 
figure  2  that  it  must  then  pass  through  two  saddle-points  with 
a  right  angle  turn  at  the  lower  one.  The  contribution  from 
the  saddle-point  at  Tg  is  subdominant  so  that  if  we  retain 
only  the  dominant  term  in  the  expansion  arising  from  then 
we  have 
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u3(x)  -  +  e4*1/6  (ic/3)1^2  x'5/12  exp(+  |  x^7^6) .  (17) 

In  an  exactly  similar  way  we  obtain 

Ug(x)  ~  -  e”711^  (‘it/3)1/2  x-5/12  exp(+  £  x^x7^) .  (18) 


Consider  next  the  solution  associated  with  the  path 
L^.  This  path  is  not  directly  deformable  into  a  steepest 
descents  path  and  we  must  therefore  choose  an  equivalent  path 
that  is  so  deformable.  A  suitable  equivalent  path  is  seen 
from  figure  2  to  be  +  B 25 *  ea°h  which  is  deformable 

into  a  steepest  descents  path.  Each  of  these  paths  must  again 
pass  through  two  saddle-points  with  a  right  angle  turn  at  the 
lower  one.  If  only  the  domaniant  term  is  retained  then  we 
have 


u4(x)  ~  +  U/3)1//2  x“5/12  exp(+  |x7/6) . 

exp(+  ^x7//^) . 


In  a  similar  way  we  obtain 

u5(x)  ~  -  (u/3)1'2  x'5//' 


(19) 

(20) 


Prom  these  results  we  can  now  define  a  set  of  standard 
solutions  that  are  real  and  unbounded  as  x  — >  +00  .  The  form 
of  the  asymptotic  expansions  (17)  to  (20)  suggests  the  follow¬ 
ing  choice : 


Bx(x)  =  ~  (u4-u5), 

b2(x)  =  2ii  (u3+U6^ 
b3(x)  =  g!  (u3-u6) . 


(21) 


and 
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The  leading  terms  In  the  asymptotic  expansions  of  these  solu¬ 
tions  are  then  given  by 


B1(x)  ~  (37t r1/2  x-5A2  exp(+  |  x7/6)# 

B2(x)  -  (3X)"1/2  x"5A2  exp(+  ^  x7//6)sln(ilZ2  x7/6  +  ^(22) 

and 

Bg(x)  ~  (3x)"ly/2  x"5/12  exp(+  2  x7//6)COB(i^2  x7/6  + 


Having  thus  tentatively  defined  a  set  of  standard  solu¬ 
tions  with  the  desired  asymptotic  behavior  as  x  — >  +oo ,  we 
turn  now  to  a  study  of  their  behavior  as  x  — >  -oo  . 
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5.  The  asymptotic  expansions  of  the  solutions  as  x  — »  -oo  . 

In  discussing  the  behavior  of  the  solutions  as  x  — >  -oo 
It  Is  convenient  to  first  replace  x  by  -x  in  equation  (4)  so 
that 

ur(-x)  -  J  exp(-xt-t7/7)dt  (s=r+l) .  (23) 

^rs 

We  again  let  t=x1^T  so  that 

ur(-x)  =  x1/6  J  exp{x7/6w(T)}dT  ,  (24) 

Lrs 

where  w(t)  now  has  the  meaning 

w (t)  «  -t-t7/7.  (25) 

The  saddle-points  therefore  still  lie  on  the  unit  circle  in 
the  x -plane  but  at  the  positions  given  by  the  six  roots  of  -1 
(see  table  2) .  The  discussion  of  the  level  curves  and  steepest 
paths  through  these  saddle  points  then  proceeds  as  before  and 
leads  to  the  qualitative  picture  shown  in  figure  3.  It  is 
useful  to  note,  however,  that  the  relative  heights  of  the 
saddle-points  are  now  given  by  the  following  scheme : 


T2 

and 

V 

+1 

T1 

and 

0 

(26) 

T 

0 

and 

V 

-1 

Consider  now  the  solution  A^-x)  which  is  still 
associated  with  the  path  L12.  11113  Path  ls  no  long©**  directly 
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deformable  Into  a  steepest  descents  path,  and  It  Is  seen  from 
figure  3  that  we  must  therefore  choose  the  equivalent  path 


L12  L15  +  L52' 


(27) 


where  the  paths  and  are  so  deformable.  They  each 
pass  through  three  saddle-points  and  make  right  angle  turns 
at  two  of  them.  The  contributions  to  the  asymptotic  expansion 
of  u^(x)  from  the  saddle-points  at  and  are  dominant 
compared  to  the  contributions  from  the  other  saddle-points . 

Thus  we  have 

y-x)  -  (3 x)-1/2  x'5712  exp(+  2^2  x7/6)sin(2  x7/6  +  (28) 


where  we  have  omitted  the  neutral  contributions  from  the 
saddle-points  at  and  and  the  subdominant  contributions 

from  tq  and  . 

Consider  next  the  solution  A g(-x)  associated  with  the 
paths  and  .  These  paths  are  directly  def ormably  into 
steepest  descents  paths,  each  of  which  then  passes  through  two 
saddle-points  with  right  angle  turns  at  one  of  them.  The 
saddle -points  at  and  therefore  make  the  same  contribution 

(except  for  a  change  of  sign)  to  Ag(-x)  as  they  do  to  A^(-x) 
and  thus  we  have 

Ag(-x)  ~  -Ax(-x).  (29) 

This  does  not  mean,  of  course,  that  A1(-x)  and  Ag(-x)  become 
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linearly  dependent*  as  x  — >  +oo ,  but  only  that  they  differ  by 
an  amount  that  is  exponentially  small  compared  to  either  of 
them.  To  investigate  this  situation  more  precisely,  consider 
their  sum  which  can  be  written  in  the  form 

A^-x)  +  Ag(-x)  =2x1  f  exp(-xt-t7/7)dt.  (30) 

L73 

Hie  major  contributions  to  ( 30)  come  from  the  saddle-points 
at  and  t^,  the  contributions  from  and  being  sub- 
dominant  by  comparison,  and  yield  a  neutral  expansion,  the 
leading  term  of  which  is 

,  A1(-x)  +  Ag( -x)  ~  (3 x)"1/2  x"5/12sin(|  x7//6  +  $).  (31) 

This  situation  does  not  occur  for  A^.C-x),  since  the  direction 
of  integration  along  one  of  the  paths  is  reversed,  and  we  have 
immediately  the  dominant  term 

A^(-x)  ~  (3x)“1//2  x“^//12exp(+  x7//^)cos(^  x7//^  +  ^) .  (32) 

This  solution  therefore  has  the  same  exponential  behavior  as 

A1(-x)  but  differs  from  it  by  n/2  in  phase. 

Turning  now  to  the  B-type  solutions  we  encounter  no 

further  difficulties.  In  the  case  of  E1(-x),  for  example,  we 

find  that  it  is  subdominant  (since  only  the  saddle-points  at 

t  and  tc  are  involved),  with  leading  term 

o  o 

— — - - - - - - 

Numerically,  of  course,  they  bocome  indistinguishable,  and 
this  is  confirmed  by  the  existing  tables  (Hughes  &  Reid  1961) . 
The  set  of  bounded  solutions  Ajc(-x)(i=l,2,3)  do  not  there¬ 
fore  form  a  "numerically  satisfactory"  set.  A  method  for 
avoiding  this  difficulty  is  discussed  in  section  7. 
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B1(-x)  ~  (3rc )-1/2  x”5/12  exp(.  x7/6)oos(^  x7//^  +  (33) 

Similarly  Bg(-x)  and  B^(-x)  are  neutral  (because  of  the  contri¬ 
butions  from  and  t^),  with  leading  terms 

B2(-x)  -  -  On)’1/2  x-5/12  sln(|  x7/6  +  f) 

and  Bg(-x)  *  +  On)-1/2  x"^/12  coa(y  x7/^  +  . 

On  comparing  these  results  with  equation  (31)  we  see  that 

B2(-x)  ~  -  A1(-x)  -  Ag(-x)  (35) 

or,,  more  precisely,  that 

A2(-x)  +  Ag(-x)  +  Bg(-x)  =  2^1  ,  exp(-xt-t7/7)dt  (36) 

L64 

and  this  is  asymptotic  to  the  subdominant  term 

(3tc)  '1/2  x-5A2  exp(-  3V3  x7A  +  5S) .  (37) 

©ius,  the  combination  of  solutions  (36)  has  the  same 
exponential  behavior  as  B^(-x)  but  differs  from  it  by  x/2  in 
phase. 
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6.  The  power-series  for  the  standard  solutions. 

In  order  to  obtain  the  power  aeries  representation  for 

the  standard  solutions  defined  in  section  4,  we  now  consider 

a  path  Lrs  (with  s  not  necessarily  equal  to  x+l)  to  consist 

of  two  rays  from  the  origin  which  bisect  the  sectors  r  and  8 . 

The  sense  of  the  path  is  from  infinity  to  the  origin  along  the 

ray  that  bisects  sector  r  and  from  the  origin  to  infinity 

along  the  ray  that  bisects  sector  s.  If  we  now  define  the 

functions  Ik(x)  by  the  relation 

oo  exp(2xik/7) 

IkU)  "  2x1  }  exp(xt-t7/7)dt  (38) 

o 

(k»0,l,2, . . .,6)  then  the  standard  solutions  can  be  expressed 
in  the  forms 

A1  =  I3"i4 

Ag  =  I2-I3+I4-I5  r  (39) 

a3  =  i(i2-i3-i4+i5) 

and 

B1  - 

b2  ‘  VW^  l  (W) 

b3  -  1<11-I2-15+i6) •  _ 

Prom  the  definition  of  Ik(x)  14  is  evident  that 

Ik(x)  *  e2*11^7  I0(xe2ltlk/7),  (4l) 

and  this  relation  can  be  used  to  express  all  of  the  functions 
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Ik(x)  (for  k^O)  appearing  in  equations  (39)  and  (40)  in  terms 
of  I0(x)  which  can,  in  turn,  he  expressed  as  the  infinite 
series 


Io(x)  -  m  7'6/7  J0{Sf]  '■  izFTs^  •  c*2) 

By  using  these  results  we  can  now  relate  the  standard  solutions 
to  the  power-series  yn(x)  defined  in  section  2  in  the  follow¬ 
ing  manner: 

6 

A-.  (x)  -  2  a  sin(6rm/7)y_(x) 

1  n«l  n  n 

6  . 

A5(x)  «  2  an{sin(4nx/7)  -  sin(6nx/7)}yn(x)  (  (43) 


n=l 

6 


n' 


U(x)  »  2  a  fcos(4nx/7)  -  cos(6rat/7)}yn(x) 
3  n=l  n  n 


and 


B-^x) 

Bg(x) 

B3W 


6 

2 

n«l 

6 

2 

n=l 

6 

2 

n-1 


an{l  -  cos(2i»t/7)}yn(x) 
an{sin(2nit/7)  -  sin(4nrc/7)}  yn(x) 
an{cos(2rw/7)  -  cos(4rut/7)}  yn(x) 


L  (44) 


'J 


where  the  coefficients  an  are  given  by 


a 


n 


7(n-7)/7 

4(n-ITr 


(45) 


It  is  also  Interesting  to  note  that  the  particular  combinations 
of  these  solutions  considered  in  equations  (30)  and  (36)  are 
given  by  the  series 
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and 


(46) 


6 

A1(x)+A2(x)  -  2  8ln(4nti/7)yn(x) 

6 

A1(x)+A2(x)+B2(x)  -  2^  an  sin(2nx/7)yn(x)  .J 

The  Wronskian  of  the  solutions  (43)  and  (44),  which  is, 
of  course,  a  constant,  is  given  by 

W ( A^ , Ag , A^,B^ , Bg , (47) 

Some  relations  which  are  useful  in  checking  the  numerical 
values  of  the  constants  appearing  in  equations  (43)  and  (44) 
are 


6 

n 

n=l 

Hi 

N 

8X3 

7  ' 

sin  sin  ^ 

sin  ^  * 

3  . 

cos  ^y  cos  ^y 

cos  = 

*  .  ' 

sin  ^y  +  sin 

V 

sin  ^  = 

!  4  • 

COS  -y  +  COS 

cos  ^  = 

1 

2 

(48) 
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7.  Discussion 

x)  and  Bk(x)  (k-1,2,3) 
defined  In  section  3  were  the  ones  tabulated  by  Hughes  &  Reid 
(1961)  and  used  by  Duty  &  Reid  (1961)  in  their  study  of  the 
stability  of  Couette  flow.  Although  the  results  obtained 
would  seem  to  be  entirely  satisfactory,  it  1b  now  clear  from 
the  discussion  given  in  section  5  (cf .  equation  (29)  in 
particular)  that  the  bounded  solutions  A^(x)  do  not  form  a 
"numerically  satisfactory"  set  of  solutions  for  large  negative 
values  of  x. 

This  situation  can  be  improved,  however,  if  we  replace 
Ag(x)  by 

I2(x)  =  A2(x)  +  A2(x) .  (49) 

The  dominant  term  in  the  asymptotic  expansion  of  Ag(x)  as 
x  — >  +ao  is,  of  course,  the  same  as  for  A2(x)  but,  as 
x  — >  -00,  Ag(x)  is  asymptotically  neutral  and,  hence,  the 
set  of  solutions  k^/K^,k^  do  form  a  numerically  satisfactory 
set.  Furthermore,  if  a  complete  set  of  six  solutions  are 
required  then  Bg(x)  is  no  longer  satisfactory  either  (cf . 
equation  (35))  and  we  should  instead  use 

l2(x)  =  X^x)  +  Bg(x) .  (50) 

The  asymmetrical  manner  in  which  the  solutions  and  Bg  have 
been  defined  here  would  appear  to  be  unavoidable  if  we  wish 
to  have  a  set  of  six  solutions  that  are  numerically 


The  standard  solutions  A.  ( 
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satisfactory  over  the  entire  Interval  -oo  <  x  <  +oo . 

Finally,  it  should  he  remarked  that  because  of  the 
oscillatory  behavior  of  these  solutions  it  would  be  advanta¬ 
geous  to  develop  the  equivalent  representation  of  them  in 
terms  of  phase  and  amplitude  functions  along  the  lines 
discussed  by  Miller  (1946)  for  the  Airy  functions. 
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Table  1.  The  saddle-point  configuration  for  x  >  0. 


Saddle-point 

Angular  argument  of 
the  position  of  the 
saddle-point 

Angle  which  the  paths 
of  steepest  descent 
make  with  the  positive 
real  T-axis 

To 

0 

0  and 

X 

T1 

x/3 

x/6  and 

-51/6 

t2 

2x/3 

it/3  and 

-2x/3 

%3 

X 

+x/2 

T4 

-2*/3 

2t/3  and 

-V3 

t5 

-it/3 

5i/6  and 

-V.6 

Table  2.  The  saddle-point  configuration  for  x  <  0. 

Saddle-point  Angular  argument  of  Angle  which  the  paths 

the  position  of  the  of  steepest  descent 
saddle-point  make  with  the  positive 


real  -c-axis 

x/6 

7it/12 

and  -5x/12 

T1 

x/2 

and  -x/4 

Tg 

5x/6 

llx/12 

and  -x/12 

t3 

-5x/6 

x/12 

and  -llx/12 

T4 

-x/2 

x/4 

and  -3x/4 

-x/6 

5x/12 

and  -7k/12 
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FIGURE  2  The  saddle-point  configuration  for  x>0 


